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Abstract 

Background and Objective: We present a new theory of selective thermal damage of 
5 non-uniformly pigmented structures in biological tissues. Spatial separation of the 
heavily pigmented areas and the target requires limitation of the pigment temperature 
and heat diffusion from the pigmented to the targeted areas. 

Study Design/Materials and Methods: A theoretical model of selective target damage 
by heat diffusion is presented for three target geometries: planar, cylindrical, and 
10 spherical. An in vitro experiment is described in which the dependence of thermal 
damage on pulsewidth at constant fluence was evaluated. 

Results: The in vitro experiment showed that the size of the damage zone for similar 
hair follicles was pulsewidth independent over a very broad range of puls'ewidths (30 
to 400 ms). We formulated a new theory (extended theory of photothermolysis) to 
15 interpret the experimental results. 

Conclusions: Based on this new theory, the treatment pulsewidth for non-uniformly 
pigmented targets is significantly longer than the target thermal relaxation time. The 
theory provides new recommendations for photoepilation and photosclerotherapy 
parameters. 

20 

Introduction 

For many years electromagnetic radiation (EMR) from lasers, lamps and other 
EMR sources (including microwave ones) have been used to treat a variety of 
medical conditions in ophthalmology, dermatology, urology, otolaryngology and other 
25 specialties. For example in dermatology EMR sources have been used to perform a 
wide variety of procedures including hair removal, treatment of various pigmented 
lesions, removal of unwanted veins, tattoo removal, and skin resurfacing. For all 
these treatments, a natural or artificial chromophore present in the body is heated by 
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absorption of either monochromatic or broadband EMR. Typical natural 
chromophores include water, melanin, hemoglobin, protein, lipid, etc.. Artificial 
chromophores can include dye, ink, carbon particles, etc. For example, heating of a 
chromophore may result directly in the destruction of a tattoo or a pigmented lesion. 
5 In these cases, the treated target for destruction and the chromophore occupy the 
same area. These cases are very well described by the theory of selective 
photothermolysis [1]. The theory of selective photothermolysis provides a quantitative 
description of optical treatments such as those mentioned above. The aim of 
selective photothermolysis (SP) is to provide permanent thermal damage of targeted 
10 structures with the surrounding tissue held intact. To satisfy the above criterion, the 
EMR pulsewidth x must be small compared to the thermal relaxation time (TRT) of 
the whole target. Actually, if the condition t«TRT is met, the heat generated within 
the target due to EMR absorption does not flow out of the structure until it becomes 
fully damaged (coagulated, injured). The above approach provides both selective 
15 damage and minimum light energy deposition. 

However, the use of such a short pulsewidth becomes inapplicable when the 
target absorption is non-uniform over its area and part of the target exhibits weak or 
no absorption, but the other part exhibits significant absorption. If this is the case, the 
weakly absorbing part of the target has to be damaged by heat diffusion from the 
20 highly pigmented/strongly absorbing one (hereafter called the heater or absorber). An 
example of such a target is the hair follicle [2]. The highly absorbing areas consist of 
melanin-bearing structures, (i.e., the hair shaft and the matrix cells). The other 
follicular tissues including the stem cells do not contain an appreciable amount of any 
chromophore that absorbs in the red/near-IR. These targets (tissues) can be 
25 damaged by heat diffusion from the hair shaft or the matrix cells to the surrounding 
follicular structures. Another example is the treatment of telangiectasia or leg veins 
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with a wavelength near the maximum hemoglobin absorption. Permanent closure of a 
vascular malformation or a vein probably requires coagulation of the vascular or vein 
wall [3]. In this case, coagulation of the wall requires heat diffusion from blood to wall. 
In this paper we will demonstrate the heat diffusion concept by presenting the 

5 results of a simple in vitro experiment involving, a hair follicle. Then, we will discuss 
the mechanism and the requirements for selective damage of biological structures by 
heat diffusion. We will prove these new concepts using thermal diffusion theory for 
isolated plane, cylindrical, and spherical targets. We will apply this theory to high- 
density targets. Finally, as an example, we will apply this new theory to the cases of 

io hair follicle and spider vein treatment. For readers who are interested in the 
mathematical details, we have included all theoretical derivations in the Appendix. 

The goal of the following experiment was to evaluate the effects of pulse 
duration and irradiance for a constant fluence on different anatomical structures of 
15 the hair follicle. The thermal damage after laser exposure was evaluated by 
histochemical demonstration of lactate dehydrogenase (LDH) activity. 

Material and methods 

Post-mortem scalp skin of one donor with light skin and dark brown hair was used for 
the experiments. The sample was procured 8 hours post mortem, frozen and stored 

20 at -80°C. The average hair density was 244 ± 26 hairs/cm 2 , and the average hair 

diameter was 68 ± 12 p.m. The hair of the post-mortem scalp sample was clipped to a 
length of approximately 0.5 mm. From this sample, 6-mm punch biopsies were 
obtained each containing an average of 69 ± 7 hairs. The biopsies were placed into 
human, hairless, post-mortem back skin of 3x3 cm 2 size and 1 .5-cm thickness that 

25 was serving as a frame. In order to obtain good fit, an opening was prepared in the 
back skin using a 5-mm skin punch. The entire sample was placed on a temperature- 

5 
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controlled plate. A 120-nm thermocouple was located below the biopsy to monitor the 
temperature of the specimen. The exposure was performed when an equilibrium 
temperature of 30±1°C was measured. 

The exposures were performed with an 800-nm diode laser prototype with 9- 

5 mm spot size and a sapphire window maintained at room temperature. The fluence 
was kept constant at 35 J/cm 2 for all exposures. The pulsewidth was 30, 100, 200, 
350, 400 and 500 ms equivalent to an irradiance of 1 170, 350,175,100, 90 and 70 
W/cm 2 . Before and after each exposure, a dermatoscopic image with crossed 
polarizing filters was taken (Dermaphot, Heine, Germany). 

10 Histochemical technique : The samples were processed in horizontal frozen 

sections of 15-u.m thickness, and a histochemical evaluation of lactate 
dehydrogenase (LDH) activity was performed as decribed by Balogh et al. [4]. 
Representative micrographs were taken at a level close to the bulge. 

15 Results 

Evaluation of the dermatoscopic images demonstrated a strong dependence 
of pulsewidth on the alterations of the clipped hair shafts. For30-ms pulsewidth, the 
hair shafts were ablated, and a brownish chromophore was released from the ablated 
hair shafts (browning). For 100-ms pulsewidth, only coarse and adjacent hair shafts 

20 demonstrated some alterations such as apparent thickening and darkening. For 
exposures with pulsewidths of 200 ms and longer, no hair-shaft alterations were 
observed at the surface after laser exposure. 

Micrographs demonstrating LDH activity at the bulge level are shown in 
Figures 1 a-e. Hair follicles with similar hair shaft diameter were selected. The inner 

25 root sheath (IRS) and the hair shaft (HS) do not possess LDH activity at the bulge 
level and therefore do not stain blue. Thermal damage of the outer root sheath (ORS) 
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or perifollicular cellular structures can be assessed by a sharply demarcated pale 
zone with loss of LDH activity. Selective loss of LDH activity of the ORS was 
observed for all pulsewidths from 30 to 400 ms. Some hair follicles also had a thin 
perfollicular cuff with loss of LDH activity. No significant variation of damage-zone 

5 diameter was observed with increasing pulsewidths for damaged single hair follicles 
at a fixed fluence of 35 J/cm 2 The extent of the thermal damage for each hair follicle 
was correlated to the coarseness and pigmentation of the hair shaft. However, for 
pulsewidths longer than 200 ms, the percentage of damaged hair follicles decreased 
with increasing pulsewidths. For a pulsewidth of 500 ms corresponding to an 

10 irradiance of 70 W/cm 2 , no decrease of LDH activity at the bulge level was observed 
for any structure. For hair follicles that were located within a few hundred microns of 
each other, inter-follicular tissue damage was observed for all pulsewidths. This 
interfollicular damage was related to hair shaft diameter, pigmentation, distance and 
local hair follicle density. 



15 



Theory and discussion 



As we can see from the experiment discussed above, the damage-zone size 
for similar hair follicles is pulsewidth independent over a very broad range of 

20 pulsewidths (30 to 400 ms). Based on the classical theory of selective 
photothermolysis [1], this behavior is very typical for pulsewidths significantly shorter 
than the target TRT. However, in our case, the pulsewidths are close to or longer 
than the TRT (not shorter!). The TRT of a hair follicle with a diameter of 200-300 jim 
is 25-55 ms; our 400-ms pulsewidth is approximately a factor of 7 longer than the 

25 maximum TRT! The next surprising result is that, even for relatively high-density 
targets such as scalp hair follicles, the thermal damage is still selective for 

7 
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pulsewidths where t»TRT. This experimental result stimulated us to investigate the 
thermal damage process of a target (such as a hair follicle) that has a structure with 
spatial separation between part of the target and the pigmented area. As we show 
below, a new theory (the extended theory of selective photothermolysis) was required 
5 to correctly interpret our experimental results. 



1. Method of selective target damage by heat diffusion 

Thermal damage of this type of target (with separation between part of the 

10 target and the pigmented area) requires deposition of sufficient heat energy into the 
absorbing area and good heat exchange between this area and the targeted external 
structures. Heat deposition depends on the absorption coefficient of the absorber and 
the EMR power density. Heat exchange depends in turn on the distance between the 
heater and the outermost part of the target and on the heat transmission coefficient 

15 between the absorber and the intervening tissue. However, at a sufficiently high 
temperature both the heater absorption coefficient and the heat transmission 
coefficient from the heater to the other targeted tissues may become lower due to 
phase transitions and destructive processes such as bleaching, melting, boiling, and 
bubble formation. This results in inefficient use of EMR energy for phase transition 

20 processes within the absorber and the intervening tissue. To prevent these 
undesirable effects, the heater peak temperature has to be limited to a prescribed 
value, Ti max , called hereafter the collapse temperature. The collapse temperature 
Timax of most natural chromophores (e.g. melanin, hemoglobin,, and water) exceeds 
100°C. Simultaneously, to ensure permanent damage of the whole target, the 

25 temperature should exceed a second prescribed value, the damage temperature T 2 , 
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throughout the target area. This temperature is lower than the collapse temperature 
Timax- More precisely, the damage temperature, T 2 , is the temperature at which 
irreversible thermal damage of the target occurs. We suggest that the basic damage 
mechanism in soft tissue is the denaturation of proteins. As an alternative to the 

5 rigorous Arrhenius rate process integral [5], the dynamic denaturation process may 
be described approximately by using the simpler damage-temperature concept. For 
human-skin collagen, T 2 is about 65-75°C and for cells T 2 is about 60°C if the 
exposure duration is several tens of milliseconds [5]. Furthermore, the tissue 
temperature should not exceed the water boiling point to prevent formation of vapor 

10 bubbles that could insulate the absorbing area from the surrounding tissue. Strictly 
speaking, vapor bubbles can transfer heat from the absorbing area to the target, but 
this process is unpredictable because vapor bubbles can move in different directions. 
To meet the temperature limitations described above, the EMR power must be 
limited, (but sufficient to heat the target up to the damage temperature) and, 

15 therefore, the EMR pulse has to be made sufficiently long to deliver enough energy. 
We define the thermal damage time (TDT) of the target to be the time required for 
irreversible target damage with sparing of the surrounding tissue. For a non-uniformly 
absorbing target structure, the TDT is the time when the outermost part of the target 
reaches the target damage temperature T 2 through heat diffusion from the heater. 

20 Because the heat-diffusion front becomes blurred when propagating into the tissue, 
part of the heat energy will leave the target. Therefore, the heated area will be larger 
than the damaged area. However, we demonstrate below that target damage can still 
be selective even though the TDT is many times as long as the TRT of the whole 
target. Apparently, the optimum EMR pulsewidth, to, should be shorter than or equal 

25 to the TDT. So, in contrast to the standard theory of selective photothermolysis, non- 
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uniformly pigmented structures have to be treated by long EMR pulses: the 
pulsewidth must typically be longer than the target TRT. In addition, the EMR power 
should be limited to prevent reduction of the heater absorption by bleaching, 
vaporization, etc. 

5 Figure 2 shows the general structure of a "target" that is thermally denatured 

by heat diffusion from a "heater". The heater includes a natural or artificial 
chromophore (pigment) with a high photon absorption coefficient. The target exhibits 
weak EMR absorption. The distance between the heater and the target is d. Photons 
from the EMR source are absorbed by the heater. As discussed in the previous 

10 section, the EMR power density has to be adjusted so that during the time of 
treatment the heater temperature T, does not exceed the temperature T 1max where 
the absorption coefficient may begin to drop. The heat propagates from the heater to 
the target due to either thermal diffusion or two other possible mechanisms. The two 
other mechanisms are hot ablation products emanating from the heater or steam. In 

15 this paper, we concentrate exclusively on the thermal diffusion mechanism. The 
process of target thermal damage is finished when its temperature reaches T z , while 
the temperature of the surrounding tissue remains below its damage temperature. In 
general, it is not possible to damage the target without damaging the tissue between 
the target and the heater. For this target type, precise selective damage of the target 

20 is impossible. But if the treated target is a layer or a shell of tissue roughly 
symmetrical with respect to the heater, the target damage can be made precisely 
selective. This means that the damage zone covers the whole target without 
extending beyond it. After the end of the EMR pulse, the target temperature at the 
" outermost point goes on growing until it reaches a maximum. The time delay 

25 between the end of the EMR pulse, t, and the moment when the peak temperature is 



JAN 25 2001 15:15 



781 993 2424 



PAGE. 11 



01/25/01 THU 15:26 FAX 781 993 2424 



PALOHAR 



12)012 



reached maximum (the latter is the TDT) is denoted by e and thus to= TDT - e. 
Therefore, the EMR pulsewidth is equal to or shorter than the TDT. This delay e is 
roughly equal to the propagation time of the heat front from the heater to the target 
e = d 2 IK where k is the thermal diffusivity of the tissue (k= 0.1 mm 2 /s); e is therefore 

5 close to the target TRT. As we shall see subsequently, TDT is significantly longer 
than the TRT in most cases. So, for selective treatment the pulsewidth should be 
approximately equal to the TDT: to=TDT. The TDT depends in turn on J^ i T 2> and 
the target size. For a better understanding of this dependence we first look more 
closely at heat diffusion from the heater to the target. Specifically, we focus on a 

io target comprising a heavily pigmented long cylinder of diameter di and a surrounding 
treated area of diameter d 2 (Fig. 3). This simple geometry can be used to model 
thermal damage of hair follicle by hair shaft heating or blood vessel destruction. We 
consider two heating modes. The first mode utilizes a "rectangular EMR pulse" (Fig. 
4a) and the second one utilizes a "rectangular temperature pulse" (Fig. 4b). In the 

is case of the rectangular EMR pulse, the heater temperature grows during the EMR 
pulse and reaches J\ at the end of the pulse (Fig. 4c). In the case of the rectangular 
temperature pulse, the temperature of the heater is constant during the EMR pulse, 
which requires a special pulse shape (Fig. 4d). For both heating modes, the heater 
temperature is below the collapse temperature Ti max so the absorption coefficient of 

20 the heater does not change. 

The sequence of thermal profiles during the heating process is depicted in 
figure 3. The input parameters for modeling are d, = 70 urn, d 2 = 210 urn. The 
collapse temperature Ti is 100°C (the boiling point of water). The damage 
temperature is 70°C, which is the average protein denaturation temperature in the 

25 10-1000 ms pulsewidth range. We assume that the EMR absorption is confined to 



I 

JAN 25 2001 15:15 



781 993 2424 PAGE. 12 



01/25/01 THU 15:26 FAX 781 993 2424 



PALOMAR 



@013 



the heater and that the thermal properties of the heater, the target and the 
surrounding tissue are the same. Heat diffusion from the heater takes place 
simultaneously with the growth of the heater temperature due to light absorption. This 
process is described well by the heat conduction equation. Figures 3a and 3b show 
5 the temperature profiles in the cylindrical target and the surrounding tissue at 
different instants of time for the rectangular EMR pulse and the rectangular 
temperature pulse, respectively. Curve 1 in both figures is the temperature profile at 
the time equal to the thermal relaxation time of the whole target. The latter time is 
TRT=d 2 2 /16k. In our case the TRT=27 ms. We can see that at the time instant t=TRT 
10 the boundary temperature of the target is still significantly below the damage 
temperature. Curve 2 in both figures is the temperature profile at the moment when 
the boundary temperature of the target reaches the damage temperature T 2 . At this 
moment, the whole target becomes damaged but the surrounding tissue is still intact. 
It is this time that has been defined above as the thermal damage time (TDT). In our 
15 case TDT=1 .8 s for the rectangular EMR pulse and TDT=360 ms for the rectangular 
temperature pulse. Based on this example, we can make two main conclusions. 
First, the ratio TDT/TRT is about 67 and 13 for the rectangular EMR pulse and the 
rectangular temperature pulse modes, respectively. Therefore, in both modes the 
pulsewidth x=TDT is significantly longer than the TRT of the entire target. Second, at 
20 the time instant when t=TDT the heated area is significantly larger than the damaged 
target. These observations present a striking contrast to the classical case of 
selective photothermolysis. The distinction above is a result of the spatial separation 
of the heavily pigmented and the treated areas. Actually, in contrast to the classical 
case, the basic damage mechanism is now heat diffusion rather than direct heating 
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by EMR absorption. The heat diffusion front is not sharp and, therefore, heat is 
spreading outside the damaged area; however, the damage is still rather selective. 

We have described the difference in treatment modality between uniformly and 
non-uniformly pigmented targets. The new extended theory of selective thermal 
5 damage of non-uniformly pigmented structures in biological tissue postulates the 
following: 

1. The EMR wavelength should be chosen to maximize contrast between the 
absorption coefficient of the pigmented area and that of the tissue surrounding the 
target. This postulate is identical to the case of. classical selective 

10 photothermolysis. 

2. The EMR power should be limited to prevent absorption loss in the pigmented 
area, but it must be sufficient to achieve a heater temperature higher than the 
target damage temperature. 

3. The pulsewidth should be made shorter than or equal to the thermal damage time 
15 (TDT), which can be significantly longer than the thermal relaxation time of the 

target. 

2. Treatment parameters for planar, cylindrical and spherical targets 

Heat diffusion is strongly dependent on the heater and target geometry. In this 
section, we discuss three basic geometries: planar, cylindrical and spherical (Fig. 5). 

20 In all cases, we assume a heater with size d! located in the center of a target with 
size d 2 . We define the ratio x= dg/d, to be the geometrical factor of the target. As 
before, Ti is the maximum pigmented area temperature, and T 2 is the target damage 
temperature (T, > T 2 ). As we show in the Appendix, the target thermal damage time 
can be expressed by the following formula: 

25 TDT=TRTxr(x,A), 
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where r(x,A) is a function of the geometrical factor x and temperature factor A defined 
as A=(T 2 - T 0 )/(Ti- T 0 ). T 0 is the target and heater temperature before irradiation. 
Normally T 0 is body temperature and is equal to 37°C. TDT is proportional to the 
TRT. In the Appendix, we present formulas for the TDT of planar, cylindrical and 
5 spherical targets. Figure 6 shows the ratio r(x,A)=TDT:TRT as a function of 
geometrical factor x for two heating modes: rectangular EMR pulse and rectangular 
temperature pulse. The calculation Parameters were Ti=100°C, T 2 =70°C, T 0 =37°C 
(A=0.52). 

We emphasize that in the framework of our analytic theory, the ratio r(x,A) 
10 does not depend on the size of the entire target and the tissue thermal properties. 
Several important conclusions follow from Figure 6. 

First, the ratio TDT/TRT is an increasing function of geometrical factor x. 
Second, the actual value of this ratio is very different for plane, cylindrical, and 
spherical targets. For a plane target, the TDT is several times higher than the TRT. 
15 The TDT exhibits appreciable growth with increasing the target dimensionality. It is 
implied herein that the planar, cylindrical, and spherical targets are one-, two-, and 
three-dimensional, respectively. 

Third, for the same TDT/TRT, the relative size of the damaged zone x is 
smallest for the spherical target. Next in this order is the cylindrical target. The plane 
20 target exhibits the largest damage area. The relative size x of the damaged zone 
around the heater decreases when increasing the target dimensionality. The latter 
two conclusions are intuitively apparent. Actually, conductive heating of a weakly 
absorbing tissue should proceed more effectively for a low-dimensional target. This 
"dimensionality" concept is a useful target parameter. It is also applicable to non- 
25 symmetrical targets. The temperature profile is sharper and better localized for the 
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spherical heater compared to the cylindrical one and it is better for the cylindrical than 
the planar heater. For the classical case of selective photothermolysis, the target 
geometry is not important because thermal damage is confined to the same area as 
the EMR absorption and direct heating. In our case, thermal damage due to heat 

5 diffusion is confined to an area that is distinct from the heater. The dependence of 
heat diffusion on heater geometry is very strong. 

Fourth, the ratio TDT/TRT depends strongly on heating mode. The rectangular 
EMR pulse mode (Fig. 4a) represents the gentlest heating mode because the heater 
temperature reaches maximum J^ at the end of the pulse (Fig. 4b). The ratio 

10 TDT/TRT is a maximum for this mode. The rectangular temperature pulse mode (Fig. 
4d) represents the most aggressive heating mode because the heater temperature 
reaches a maximum just after the beginning of the pulse and the maximum heater 
temperature takes place during the EMR pulse. The ratio TDT/TRT is a minimum for 
the rectangular temperature pulse mode. As mentioned above the rectangular 

15 temperature pulse mode can be realized by using an EMR pulse with a special 
temporal profile. The initial power density should be very high to raise the heater 
temperature abruptly (for a time interval on the order of or shorter than the TRT). 
After the maximum temperature is reached, the power density should undergo a 
steep fall to prevent overheating. Then, to maintain the heater temperature at the 

20 prescribed level (Fig. 4c), the power density should fade gradually to compensate the 
heat flow out of the heater. The pulse power should be precisely adjusted to keep the 
heater temperature below the collapse temperature Ti. The power depends on heater 
absorption and size and the EMR attenuation in tissue (see Appendix B). In reality, it 
is probably difficult to exactly create these two modes. Thus, the real value of TDT 

25 can be between these two extreme cases. 
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The ratio TDT/TRT depends on the temperature factor A=(T 2 - T 0 )/(T r T 0 ). 
Table 1 shows the influence of initial target temperature T 0 and maximum heater 
temperature T v All calculations were done for a cylindrical target with the same size 
as figure 4 for the rectangular EMR pulse. The TDT increases by a factor of 2.2 by 

5 precooling from 37°C to 27°C and decreases by a factor of 3 by preheating to 45°C. 
The fluence should be changed at the same time. If the heater temperature can 
reach a high value without losing absorption, the TDT can be significantly reduced. 
The ratio TDT/TRT is close to two (TDT=2TRT) for the case when the heater 
temperature is 200-250°C. In biological tissue, such a high temperature can be 

10 expected in melanin in the hair shaft or in an artificial chromophore such as carbon. 
But we must remember that the thermal diffusivity can drop in the tissue surrounding 
the heater due to water vaporization. So this case is very difficult to predict. 



15 



Table 1. Thermal damage time (TDT) as function of temperature factor 



Initial 

temperature 
To/C 


Temperature 
of heater 
Ti,°C 


Temperature 
of damage 
T 2 ,°C 


Temperature 

factor 

A 


TDT/TRT 


37 


100 


70 


0.53 


67 


27(precooling) 


100 


70 


0.59 


150 


45(preheating) 


100 


70 


0.45 


22 


37 


200 


70 


0.21 


2.5 


37 


240 


70 


0.17 


1.8 



different target geometries. The heater temperature should depend on the heater 
geometry. Figure 7 shows the heater center temperature as a function of pulsewidth 
for a rectangular EMR pulse with the same power. Spherical, cylindrical, and planar 
20 heaters have equal: size d tl thermal properties, and EMR absorption coefficient. The 
thermal relaxation time of the heater x r depends on geometry, and the ratio is 1:2:3 
for spherical, cylindrical, and planar heaters, respectively. If the pulsewidth % is 
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significantly shorter than the thermal relaxation time of the heater x, (t«x r ) the 
temperature rises of all the heaters exhibit the same elevation of temperature. 
However, as shown in Figure 6, the temperature behavior of heaters with different 
geometries is very different for pulsewidth equal to or longer than x r . A steady state 

5 heater temperature for a rectangular EMR pulse is possible only for a spherical 
heater. For cylindrical and especially for planar heaters, the temperature is 
continuously rising when the pulsewidth is increasing (the power density should be 
constant, the energy density should be proportional to the pulsewidth). This is 
because 3D heat diffusion from the spherical heater (in contrast with 2D and 1D heat 

10 diffusion from the cylindrical and planar heaters respectively) can compensate 
constant heating from the rectangular EMR pulse. This phenomenon is very 
important for uniform targets when the heater and target are the same. For example, 
the temperature of epidermis (planar target) for a rectangular optical pulse at a 
wavelength strongly absorbed by melanin will continuously rise during a long pulse (x 

15 »x r ). However, for a spherical target such as the hair bulb matrix, the temperature 
will stabilize at a steady-state level, To produce a constant heater temperature 
(rectangular temperature pulse as shown in Figure 4d, the EMR pulse shape must be 
special (Fig. 4c) with a strong peak in the beginning and decaying amplitude. The 
ratio of peak and decaying amplitudes are different for heaters with planar, cylinder 

20 and spherical geometries (Fig. 8) 

3. Mechanism of selective damage of high density targets by heat diffusion. 
Since the results obtained in the above sections were presented just for 
isolated targets, it is important to see the effect of high-density targets. As shown 
above, non-pigmented targets located close to pigmented areas can be destroyed by 
25 heat diffusion. However, heat flux will leak out of the target and the heating area will 
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be larger than the target. This increases the risk of overheating the tissue 
surrounding the target and thus the risk of unselective damage. Very roughly the 
fluence to produce unselective bulk tissue damage F NS is (d3/d 2 ) n times higher than 
the fluence required to produce selective damage Fs, where 62 and d3 are the target 
5 size and distance between target centers, and n is 1 ,2 or 3 for planar, cylindrical, and 
spherical targets, respectively. So as a first approximation, Fns /Fs is proportional to 
the ratio of target volume and tissue bulk volume and is independent of pulsewidth. 
The risk of unselective damage increases in the following order: spherical, cylindrical 
and planar targets. 

10 More precisely, consider this problem for the ideal model of cylindrical targets 

with equal spacing. Figures 9 shows the dependence of TDT, F s , and F N $ on density 
fectory=d3/d 2 for a rectangular temperature pulse. As seen from Figure 9, F s 
decreases in y beginning with y=4. This effect is explained by the influence of heat 
fluxes from neighboring targets. But at the same time, the Fns is continuously 

15 dropping as function of density factor y. So for y<4 the range of safe fluences is going 
to be very narrow and the risk of unselective damage will increase. In the case of 
short pulse/high-power irradiation where the chromophore is ablated, the energy 
delivered to the tissue is therefore limited. This provides additional protection in the 
case of high-density target treatment. 

20 4. Applications of the extended theory of selective photothermolysis 

Now we are going to discuss applications of the new theory. Thus far, we have 
considered at some length the method of selective photodamage of unpigmented 
targets in the tissue. We will present two examples where the new theory can be 
applied; photoepilation and photosclerotherapy. 
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Photoepilation. Photoepilation utilizes light to cause thermal or mechanical 
damage of hair follicles. To achieve hair growth delay, it is sufficient to either 
damage matrix cells of anagen hair follicles or coagulate blood vessels of the papilla 
or possibly to destroy part of the outer root sheath [2]. For permanent hair-follicle 
5 damage, in accordance with current knowledge, it is necessary to damage stem cells 
that are located in the bulge area at the interface of the outer root sheath and the 
connective tissue sheath [6]. One can also irreversibly damage a hair follicle at the 
level of the dermis by replacing it with connective tissue. 

The matrix cells produce the hair shaft. The matrix cells contain melanosomes 

10 that produce hair melanin. The concentration of melanin in the matrix cells is 

significantly higher than in the hair shaft. Melanin is distributed uniformly and densely 
in the matrix cells. So for a pulsewidth longer than the TRT of individual 
melanosomes (1|is), the matrix cells act as a uniformly pigmented target. This is a 
typical example of a target where standard selective photothermolysis theory is 

15 applicable. For selective and effective treatment, the energy and pulsewidth have to 
be significantly shorter than the TRT of the matrix cells. The matrix cells form a dome 
shaped structure with the smallest cells close to half of the hair shaft diameter d h . So 
the TRT of the hair matrix can be estimated as the TRT of a layer with thickness dh/2, 
that is d h 2 /32k. The pulsewidth values for effective treatment are presented in Table 

20 2. 
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Table 2. Optimum pulsewidth for hair-follicle treatment. 



Hair 

type/diameter 


TRT of 

hair 

shaft, 

ms 


TRT of 

hair 

follicle, 

ms 


Halting hair shaft growth 


Permanent hair-follicle 
destruction 


TRT of 
hair 
matrix, 
ms 


TDT of papilla blood vessel, 
ms 


TDT of stem cell, ms 


Rectangular 
temperature 
pulse 


Rectangular 
light pulse 


Rectangular 
temperature 
pulse 


Rectangular 
light pulse 


Fine/30u 


0.6 


5.4 


<0.3 


0.5 


1.5 


30 


115 


Medium 
coarse/70ji 


3 


27 


<2 


2.7 


8.5 


170 


610 


Large 

coarse/1 20m. 


9.6 


87 


<5 


8 


21 


510 


1800 



5 Another method of halting hair shaft growth is to coagulate blood vessels in 



the papilla. The loop of blood vessels in the papilla is located in the center of the 
matrix cell dome. Blood absorption is significantly lower than melanin absorption in 
the neighboring matrix cells because of the small vessel size. So the most effective 
method of papilla blood vessel coagulation is to utilize heat diffusion from the matrix 

10 cells that absorb light. This is a typical case where the extended theory of selective 
photothermolysis can be applied. The highly pigmented heater (the matrix cells) and 
the lightly pigmented target (the blood vessels) are separated by distance about d h /2. 
We calculated the TDT of the hair papilla blood vessels using our thermal diffusion 
model assuming the maximum heater temperature to be Ti=100°C and the blood 

15 vessel coagulation temperature to be T 2 =65°C. The TDT value for different hair sizes 
and heating modes are presented in Table 1 . We can see a significant difference in 
optimum treatment pulsewidth for: 1) the matrix cells by direct light absorption and 2) 
the papilla blood vessel loop by heat diffusion. 

Let us now consider damage of hair stem cells. They are located in the basal 

20 cell layer of the outer root sheath of the lower isthmus. The stem cells do not have 
any pigment that can effectively absorb light in the therapeutic window (600-1200 
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nm), which is also the best wavelength range for photoepilation [2]. However, the 
stem cells can be damaged by heat diffusion from the melanin-rich hair shaft or an 
artificial chromophore inside the inner root sheath. This case is very well described 
by the cylindrical model that we considered in some detail above. The optimum 

5 treatment pulsewidth of the stem cells is close to the TDT of the follicle structure. If , 
we assume that the maximum hair-shaft temperature is Ti=100°C and the damage 
temperature of the stem cells is Ti=65°C f we can calculate the TDT as shown in 
Table B (approximate formulas give us slightly lower values than the exact modeling 
we presented in Figure 6b). Calculation results for a follicle with a ratio of follicle 

10 diameter to hair-shaft diameter of x-3 are presented in Table 2. 

To estimate power density P and fluence F for the rectangular light pulse case, 
we can use formulas from Table 2A (See Appendix ) that we can simplify for the 
follicle case as: 



F = P'TDT[Jtcm*] 

where ^ a [cm 1 ] is the hair-shaft absorption coefficient, q is the ratio of the radiance at 
20 the target location to the input power density (attenuation factor). For our calculation, 
we use typical treatment conditions: wavelength is 800 nm, dark hair with ^ a =100 cm" 
\ x=3, q=1(low isthmus level of follicle), T l= 100 o C , T 2 =65°C, and T 0 =36°C. Using 
the formulas shown above, we calculate the following values: power density is 560 




2 



P = 
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W/cm 2 for 30 \im fine hair, 100 W/cm 2 for 70 \m\ terminal hair and 35 W/cm 2 for 
120 \im coarse hair. For these hairs with equal melanin concentration and different 
diameters, the fluence necessary to damage the stem cells is the same. The fluence 
value is 40 J/cm 2 . The pulsewidth value appears in Table 2 as the TDT for a 

5 rectangular light pulse. We can see from Table 2 that selective and complete hair- 
follicle damage can be achieved over a very broad range of pulsewidths. For 
example, for a hair shaft diameter of approximately di=70 |xm and a hair follicle 
diameter of di=210 nm, this range is 170-610 ms, which is significantly longer than 
the TRT (27 ms). In our in vitro experiment, we observed nearly complete damage of 

10 similar size follicles for pulsewidths in the 30-400 ms range with an identical fluence 
of 35 J/cm 2 . The main differences between theory and experiment for short pulses 
can be explained by differences in the maximum hair-shaft temperature used in the 
theoretical model (Tj =100°C) and the hair-shaft temperature reached in the 
experiment. For a pulsewidth shorter than 100 ms, we saw a significant change in the 

is hair-shaft structure. This implies that the hair-shaft temperature was higher than 
250°C. For this temperature the TDT is very close to the TRT (see Table 1 ). If the 
energy of the ablation products is not too high, most of them release their energy 
inside the hair follicle and heat up the follicle structure. So all the absorbed energy 
was utilized to heat the follicle structure (including IRS and ORS). For a pulsewidth 

20 significantly shorter than 30 ms with the same fluence, we can expect more hair-shaft 
ablation" and escape of the hair-shaft ablation products from the hair follicle. In this 
case, the damaged volume of the hair follicle should decrease with decreasing 
pulsewidth and, for very short pulses, it should be limited to the hair shaft. The 
differences between theory and experiment for pulses longer than 400 ms can be 
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explained by the low power density (<100 W/cm 2 ) of the laser pulse, which was not 
sufficient to heat the hair shaft close to 1 00°C. 

Photosclerotherapy. Photosclerotherapy produces thermal or mechanical 

5 damage of the vessel structure due to light absorption by blood [6]. The blood 
hemoglobin exhibits selective light absorption over a wide wavelength range. 
Optimum vessel closure can be achieved by denaturation of the endothelium that is 
in direct contact with blood. This case is well described by the standard theory of 
selective photothermolysisp]. Other authors have suggested that permanent vessel 

io closure requires denaturation of the vessel wall structure including the tunica intima, 
tunica medula, and tunica adventitial]. These structures do not contain any strongly 
absorbing chromophores and can be damaged by heat diffusion from blood. 
Therefore, we will apply the extended theory of selective photothermolysis to 
estimate treatment parameters for different vein sizes. As a first approximation, the 

15 vein can be modeled as by a cylinder. However, this is true just for limited 

constrained cases depending on the light wavelength and vein size. The cylindrical 
model is valid if the light penetration depth in the blood exceeds the vessel internal 
diameter D. In this case, the blood is heated uniformly. If scattering of the blood is 
lower than absorption, the light penetration depth is roughly equal to the inverse 

20 blood absorption coefficient 1/jx a - Thus the cylindrical model can be used when 
D<1/p* . In this case, the heater diameter is di=D and the target diameter is d 2 - 
=D+2w, where w is the vessel-wall thickness. If D»1/^a, the heated zone is a 
cylindrical blood layer, in contact with the vessel wall. The heated-zone thickness is 
approximately equal tot the penetration depth 1/|x a . This case can be roughly 

25 described by the planar model with heater thickness di=1/Ha and target thickness d 2 - 
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=1/[ia+2w. Analytic evaluation of the TDT using the expressions of Table 2A is 
possible in the rectangular light-pulse case when the whole blood volume is uniformly 
irradiated. If this is the case the cylindrical target model is applicable. In the opposite 
limiting case of very strong blood absorption, a thin layer of blood adjacent to the 

5 vessel wall is only irradiated. In the latter case one may apply the planar model. At 
that, the analytic and numeric results show rather good agreement. In other cases 
both the TDT and the input flux were evaluated numerically. The maximum heater 
temperature is T 1 =100°C (limited by blood coagulation and vaporization), vessel wall 
denaturation temperature is T 2 =:65°C, and the initial body temperature is T 0 =37°C. 

10 The calculations were performed for two wavelengths X: X=577 nm (maximum 
hemoglobin absorption), where ^/\l a -43 \xm and moderate hemoglobin absorption 
?t=1060 nm, where 1/^=1400 p,m [8,9]. To estimate the power density and fluence 
for a rectangular light pulse, we can use formulas from Table B and direct modeling 
for a rectangular temperature pulse. Calculation results for different types of spider 

15 veins are presented in Table 3. 



Table 3. Optimum pulsewidth for treatment of spider veins. 



Diameter 
of vein, 
mm 


Wall 

thickness, 
mm 


TRT of 
blood 
volume, 
ms 


TRT 
of 

vein, 
ms 


TDT of vein and fluence F 


a=577 nm 


X-1060 nm* 


Rectangular 
temperature 
pulse 


Rectangular 
light pulse 


Rectangular 
temperature 
pulse 


Rectangular 
light pulse 


TDT, 
ms 


F, 

J/cm 2 


TDT, 

ms 


J/cm 2 


TDT, 
ms 


J/cm 2 


TDT, 
ms 


J/cm 2 


0.1 


0.035 


0.5 


5 


40 


10 


130 


20 


40 


590 


140 


1090 


0.25 


0.08 


4 


35 


150 


6 


515 


15 


150 


335 


610 


615 


0.5 


0.12 


35 


135 


215 


4 


740 


6 


215 


115 


1200 


200 


1 


0.15 


260 


540 


240 


4 


670 


5 


240 


50 


2400 


90 



JRN 25 2001 15^23 



781 993 2424 PAGE. 08 



01/25/01 THU 15:34 FAX 781 993 2424 



PALOMAR 



@009 



* We can predict that TDT and F will be on the same order for another important 
wavelength (?l=800 nm) because the absorption coefficients for vein blood at 
X^=800 nm and \=1060 nm are similar. 

As we see from Table 3, the TDT of the entire vein wall structure is shorter for 

5 the 577-nm wavelength that coincides with the strong hemoglobin absorption peak. 
The TDT can be very long, typically for large veins. The perfusion factor can also be 
important but is neglected in the present study. Blocking of the blood flow can be 
important for optimum treatment with such a super long pulse. The large veins need 
signif icantlt lower fluence, power density, and longer pulsewidth than small veins. 

10 This dependence is even stronger for light with low absorption in hemoglobin. The 
water absorption of surrounding tissue can be important for long treatment 
wavelengths (^>750nm) due to parallel non-selective heating by water absorption 
that can increase the tissue temperature T 0 . As we showed above, the TDT will be 
shorter in this case and the fluence will be lower. 

is The parameters for spider-vein treatment suggested by our new theory are 

very different compared to typical clinical parameters used for photosclerotherapy 
and should be clinically proven. 

Following from tables 2 and 3, the treatment time (pulsewidth) for hair follicles 
and spider veins is on the order of several hundreds of milliseconds. As was showed 

20 in [10], parallel cooling of the epidermis (simultaneous heating by light absorption and 
heat removal by heat diffusion into the skin and cooling agent) is very effective in this 
pulsewidth range. This cooling mode is important for epidermal protection and makes 
the procedure more effective because high fluence can be delivered through the 
epidermis. 
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Conclusions 

1. Our in vitro experiment demonstrated that the hair-follicle damage pattern for fixed 
5 fluence is selective and pulsewidth independent over a broad pulsewidth range; 

this range is 1-10 times the hair-follicle TRT. 

2. A new theory of selective thermal damage of non-uniformly pigmented structures 
of biological tissues is presented in this paper. We call this theory the extended 
theory of selective photothermolysis (ESP). 

10 3. The extended theory of selective photothermolysis is important for biological 
targets with space between all or part of the target and the pigmented area of the 
tissue. So, in contrast with the standard theory of selective photothermolysis the 
target is destroyed by heat diffusion from the pigmented area to the target but not 
by direct heating from electromagnetic radiation (EMR). 

15 4. The main principles of the extended theory of selective photothermolysis are: 

• The EMR wavelength should be chosen to maximize contrast between the 
absorption coefficient of the pigmented area and that of the tissue surrounding the 
target. This postulate is identical to the case of classical selective 
photothermolysis. 

20 • The EMR power should be limited to prevent absorption loss in the pigmented 
area, but it must be sufficient to achieve a temperature of the pigmented area 
higher than the target damage temperature. 

• The pulsewidth should be made shorter than or equal to the thermal damage time 
(TDT). The TDT is the time required for irreversible target damage with sparing of 
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the surrounding tissue. The TDT can be significantly longer than the thermal 
relaxation time (TRT) of the entire target. 

5. The power and shape of the EMR pulse are important for successful treatment in 
contrast with the standard theory of selective photothermolysis where just fluence 

5 is important. 

6. Treatment parameters are power and pulsewidth dependent and are not just 
dependent on target size, absorption and thermal properties but their geometry 
as well. 

7. For highly dense targets, it is important to limit fluence to prevent non-selective 
10 damage. Selective treatment of very high density targets is not possible because 

the gap between the fluence of selective and non-selective treatment is too small. 

8. The extended theory of selective photothermolysis suggests a significantly longer 
treatment time (pulsewidth) than classical selective photothermolysis. In 
dermatology, it provides the additional benefit of epidermal protection due to 

15 parallel cooling. 

9. The extended theory of selective photothermolysis presented in this paper is 
useful for electrothermolysis (electrode is heater) and certain cases of magnetic 
and acoustic treatment of tissue. 

20 Appendix. 

Procedure to determine the amplitude and the duration of the rectangular EMR 
pulse: brief summary. 

The present Appendix summarizes the EMR pulse parameters for the 
treatment of the basic targets exhibiting a high degree of symmetry, i.e., the planar, 
25 cylindrical, and spherical ones. 
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The notations and the basic parameters of the problem are explained in Table 

1A: 



Table A1 



Variable 


Dimensionality 


Name 


Assumptions and relations 


k 


5 — T" 

cm s 


Thermal diffusivity 


Assumed to be the same all 
over the tarqet 


P 


gcm a 


Density 


Assumed to be the same all 
over the tarqet 


L> 




Specific heat 


Assumed to be the same all 
over the tarqet 


Ho 


cm" 1 


Tissue absorption coefficient 


Assumed to be zero outside the 
heater 


q 


a.u. 


The ratio of radiance to the 
input power density 




di 


Cm 


Thickness or diameter of the 
heater 




d 2 


Cm 


Thickness or diameter of the 
tarqet 


d 2 >d i 


a 3 


rim 


Mean spacing between the 
tarqets 


d > d 


To 


u c 


Initial temperature of both the 
target and the surrounding 
tissue 


r 0 = 37°C 

0 


Tlmax 


U C 


Temperature of heater 
absorption toss, collapse 
temperature 


r lira ,=ioo 0 c-25o 0 c 


T, 


°C 


Maximum temperature of the 
heater(absorber) 


T2<T X <r lmax 


T 2 


°C 


Temperature of irreversible 
damage of the tissue 


r 2 =70°c 


A 


a.u. 


Temperature factor, 
temperature ratio 


A = r2 " r ° <1 


X 


a.u. 


Geometrical factor, diameter 
ratio 


xadjdi >1 



Our analysis was based on the heat conduction equation. We have found 



approximate analytic solutions of thermal damage time (TDT) and input power 
density P. 

The final expressions for the important variables in question are outlined in 
simplified form in Table 2A below. It is implied that the thermal constants, i.e., the 
10 density, the thermal diffusivity, and the thermal conductivity do not vary significantly 
within the target and the surrounding tissues. 
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Table 2A 



N 



Quantity 



Notation 

[Dimensi 

onality] 



Approximate expression for a particular target geometry: 



Planar 



Cylindrical 



Spherical 



Thermal 
relaxation 
time of the 
heater 



x, 
[s] 



r 8-ft 



T. - 



' 16 k 



41 
' 24 k 



Thermal 
relaxation 
lime of the 
target 



TRT 
[s] 



TRT = -*- = 



TRT = 



ill 
16Jfc 



TRT =- 



■ = X ■ T, 



= X-T r 



Thermal 
damage 
time 



TDT 

[s] 



TDT = 



TRT 



TDT = 



TRT 



TDT = 



2-x 4 



D-A 
1-A 



-1 



•exp 



D-0.3 A 



\ 



D - exp(- *^)+ 
+ 1.8- x- erf (x). 



0 = 0.6 + 
+ 2-ln(x)- 
-Ei(-1.4x 2 ) 



0.9- 



TRT 



1-A 



D-A 



D-A>0 
o°, D-AZO. 



D = 0J 



erf(l.3;c) 



Input power 
density 



[W/cm 2 ] 



P = 



pc l.lV 



Mai 



TRT 



T r T o 



1 



2 TDT 

1+2.1-jT 1 

TRT 



P = 



pc 



M a 1 TRT 



P = 



p c 0.3-.X 



Hq TRT 



In 



1+1 .4-jr 



V 



2 TDT ) 
TRT 



1- 



1 



1+1. 2- x 



2 TDT 

TRT 



Input 
fluence 



[J/cm 2 1 



F = P-TDT 



F = PTDT 



F = P- TDT 



10 



The present discussion is restricted to rectangular EMR pulses only. Our goal 
herewith is to determine the pulsewidth and the pulse amplitude. This may be 
performed in the following order: 

1. Based on the target geometry, the size of the pigmented area (heater) d 1f size of 
theitarget d 2 , distance between targets d 3> the geometrical factor x= d 2 /di and 
density factor can be determined. 
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2. Based on the thermal properties of the target and chromophore, the temperature 
of pigmented area absorption loss Ti max , temperature of target damage T 2 , initial 
tissue temperature T 0 , the temperature factor A=(T 2 - T 0 ) /(TV T 0 ) can be 
determined. 

5 3. Using these parameters and formulas, we can use the formulas in Table 2A, rows 
1-3 to determine x r , TRT and finally TDT. As explained above, the TDT is 
approximately equal to the optimum duration of the EMR pulse t 0 . 

4. Based upon an estimate of the EMR attenuation factor at the depth of target q 
and the absorption coefficient of pigmented area m and the heater temperature 

10 T\, the power density on the skin P using can be determined by using the 

formulas from Table 2A row 4. The power density should be limited in order to not 
bleach the pigmented area, but it should be significantly high to rich the 
temperature of target damage T 2 . 

5. Treatment fluence is given by F = P TDTor F = P r 0 

15 

In contrast to both the planar and cylindrical targets, the TDT for the spherical 
target may be evaluated to infinity (see Table 2A, row 3, rightmost column). This 
means that one cannot ensure the safety of the heater in an attempt to damage the 
whole target. After a part of the target becomes damaged, the heater temperature 
20 reaches the crucial value T t . This gives rise to phase transitions, bleaching, bubble 
formatijon, and other nonlinear processes lying outside the scope of this paper. 
Therefore, our simple theory provides the means to describe thermal damage of 
sufficiently small spherical targets only. More precisely, for a given value of A the 
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diameter ratio x must not exceed the value obtained from the equation D = 0, where 
variable D is a function of x determined by the last expression of Table 2A, row 3, 
rightmost column. 

5 

Figures. 

Figure 1. Thermal damage pattern of hair follicle (HF) at bulge level demonstrated by 

10 loss of LDH activity, Fluence was 35 J/cm 2 . Bar, 100 urn There is selective damage 
of the hair follicle over a wide range of pulsewidths, and the extent of thermal 
damage does depend on hair-shaft diameter and pigmentation, (a) x= 30 ms, l= 1 170 
W/cm 2 . Right HF demonstrates almost complete damage of outer root sheath (ORS), 
hair shaft (HS) structurally disintegrated with browning, fragments of HS partly 

15 protruded through ruptured IRS, no perifollicular damage. Left HF presents vacuoles 
within HS, no browning, IRS ruptured, partial loss of LDH activity of ORS and no 
perifollicular damage. Middle HF with fine and light pigmented HS does not 
demonstrate any HS alteration or loss of LDH activity, (b) x= 100 ms, l= 350 W/cm 2 . 
Right HF without any HS alteration or loss of LDH activity of ORS. Left HF with 

20 vacuoles within HS, IRS ruptured, loss of LDH stain of entire ORS and a thin cuff of 
perifollicular tissue, (c) x= 200 ms, 1=175 W/cm 2 . No hair shaft alterations for both HF. 
Right HF with minor loss of LDH activity of ORS. Left HF with partial loss of LDH 
activity of ORS. No perifollicular damage, (d) z= 200 ms, 1=175 W/cm 2 . No hair shaft 
alterations for both HF. Right HF with partial loss of LDH activity of ORS. Left HF with 

25 complete loss of LDH activity of ORS. Interfollicular damage of these adjacent HF 
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with only 300 um interfollicular distance. Note more pronounced pigmentation of hair 
shafts than in fig. 1c. (e) x= 350 ms, 1=100 W/cm 2 . No hair shaft alterations for both 
HF, IRS not damaged, partial loss of LDH activity for both HF, more pronounced loss 
for HF with more pigmented HS. 

5 

Figure 2. General structure of biological target with space separation between target 
and absorber (heater). 

Figure 3. Temperature distribution in the tissue with cylindrical absorber with 
10 diameter 0.07 mm and target with diameter 0.21 mm. Figure 3a shows the 
temperature distribution for a rectangular EMR pulse at two points in time: bottom 
curve at t=TRT=27.5ms and top curve at t=TDT=1.6s. Figure 3b shows the 
temperature distribution for a rectangular temperature pulse at two points in time: 
bottom curve at t=TRT=27.5ms and top curve at t=TDT=0.36s. Maximum 
15 temperature of the absorber is Ti =100°C. Damage temperature is T 2 =70°C. Initial 
temperature is T 0 =37°C. The absorption of surrounding heater tissue was neglected. 

Figure 4. Time dependence of electromagnetic radiation (EMR) power and 
temperature of the heater (absorber). Two basic cases are shown: (Fig. 4a) 
20 rectangular EMR pulse and (Fig. 4c) rectangular temperature pulse. Figure 4b shows 
heater temperature as a function of time that is produced by a rectangular EMR 
pulse. Figure 4d shows EMR power as function of time that produces a rectangular 
temperature pulse at the heater. 
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Figure 5. Three types of target with different geometry: (a) planar, (b) cylindrical, (c) 
spherical. 1 is the heater (absorber). 2 is the target, d, is size of the heater, d 2 is size 
of the target. 

5 Figure 6. The ratio r of thermal damage time (TDT) and thermal relaxation time (TRT) 
as function of ratio x of size of the target d 2 and size of the heater di for different 
target geometries: (Figure 6a) planar, (Figure 6b) cylindrical, (Figure 6c) spherical. 
Ratio ri=TDT:TRT is given for two heating modes: rectangular EMR pulse and 
rectangular temperature pulse. 

10 

Figure 7. Increase of heater temperature T(t) normalized to steady state and 
temperature of spherical heater increase T 0 as a function of heating time for 
rectangular pulse and different type of heaters with the same size: (1) planar, (2) 
cylinder and (3) sphere. 

15 

Figure 8. EMR power P(t) pulses producing rectangular temperature pulse on the 
heater normalized to peak power P 0 at t=0 for different target geometries with the 
same size: (1) planar, (2) cylinder and (3) sphere. 

20 Figure 9. Thermal damage time and damage fluence for high density targets as a 
function of density factor y (ratio of distance between centers of targets d 3 and size of 
the target d 2 ). Figure 9a shows the dependence of TDT for selective (bottom curve) 
and nonselective (top curve) damage on density factor y. Figure 9b shows 
dependence of damage fluence for selective F s (bottom curve) and nonselective F NS 

25 (top curve) damage on density factor y. Both curves are normalized to the damage 
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fluence of individual target F 0 . Calculations were done for a cylindrical individual 
target with parameters the same as those used for Figure 3. 

Figure A1 . Heat production as a function of radius. The approximation of the stepwise 
function distribution (solid curve) by the Gaussian distribution (dotted curve). 

5 
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Figure 3b. 
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Figure 4a. 
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Figure 4c. 
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Figure 4d. 
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Figure 5 a. 



JAN 25 2001 15:30 



781 333 2424 



PAGE. 11 



01/25/01 THU 15:41 FAX 781 993 2424 



PALOMAR 



@)012 




01/25/01 THU 15:41 FAX 781 993 2424 



PALOMAR 



(21013 



di 




Figure 5c. 
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Figure 6a. 



JPN 25 2001 15=30 



781 993 2424 



PAGE. 14 



01/25/01 THU 15:41 FAX 781 993 2424 



PALOMAR 



©015 




1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 3.4 3.6 3.8 4 



x=dyd t 



Figure 6b. 
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Figure 6c. 
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Figure 7. 
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Figure 9a. 
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Figure 9b. 
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